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ABSTRACT 


Aims. We present the first detailed analysis of the astrophysical parameters of the poorly studied Sco-Cen member HD 152384 and its 
circumstellar environment. 

Methods. We analyse newly obtained optical-near-IR XSHOOTER spectra, as well as archival TESS data, of HD 152384. In addition, 
we use literature photometric data to construct a detailed spectral energy distribution (SED) of the star. 

Results. The photospheric absorption lines in the spectrum of HD 152384 are characteristic of a AO V star, for which we derive a 
stellar mass of 2.1 + 0.1 Mo and a stellar age > 4.5 Myr. Superimposed on the photospheric absorption, the optical spectrum also 
displays double-peaked emission lines of Can, Fe1, Mgt and Sit, typical of circumstellar disks. Notably, all Hydrogen and Helium 
lines appear strictly in absorption. A toy model shows that the observed emission line profiles can be reproduced by emission from a 
compact (radius < 0.3 au) disk seen at an inclination of ~ 24°. Further evidence for the presence of circumstellar material comes from 
the detection of a moderate infrared excess in the SED, similar to those found in extreme debris disk systems. 

Conclusions. We conclude that HD 152384 is surrounded by a tenuous circumstellar disk which, although rich in refractory elements, 
is highly depleted of volatile elements. To the best of our knowledge such a disk is unique within the group of young stars. However, 
it is reminiscent of the disks seen in some white dwarfs, which have been attributed to the disruption of rocky planets. We suggest 
that the disk around HD 152384 may have a similar origin and may be due to collisions in a newly formed planetary system. 
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œœ 1. Introduction 
N 
CNA HD 152384 (HIP 82714, G = 7.01) is a poorly studied early-type 
OO star located in the constellation of Scorpius 
© 2012). From its parallax and proper motions, 
Ya (2011) and [Galli et al.| (2018) identified HD 152384 
© as a member of the 5-10 Myr Sco-OB2 association, suggest- 
w= ing a young age. Based on an anomaly between the Hip- 
CN parcos and Gaia DR2 proper motions, 
XS classified HD 152384 as a suspected binary system. More re- 
«re cently, the multiple nature of HD 152384 was confirmed by 
(2021), who identified the source Gaia EDR3 
= 60262887 13470340608 (G = 13.58, separation 77.20” or 9908 
au) as its companion. 

HD 152384 is included in a new sample of 246 Gaia- 
selected young intermediate-mass stars for which we are cur- 
rently analysing UV-near-IR spectroscopy obtained with the 
XSHOOTER instrument on ESO’s Very Large Telescope (VLT). 
Within this sample, and to our knowledge within the entire class 
of YSOs, HD 152384 shows a unique optical-near-IR spectrum, 
on which we report in this Letter. 


* Based on observations collected at the European Southern Obser- 
vatory under ESO programme 0101.C-0902(A). 


2. Observations 


A high-quality spectrum of HD 152384 was obtained with 
XSHOOTER - a wide spectral range (300-2400 nm) spectro- 
graph on the VLT —on July 2, 2018. Total 
integration times were 60 seconds, 60 seconds and 100 seconds 
for XSHOOTER’s UV, visual and near-IR arms, respectively. 
The XSHOOTER entrance slits were oriented along the paral- 
lactic angle (—93.9°) to minimise atmospheric diffraction losses 
and slit widths were set to 1.0", 0.9" and 0.6" for the UV, vi- 
sual and near-IR arms, resulting in spectra with a spectral reso- 
lution of 5400, 8900 and 8100 for these three arms. Data were 
reduced using the XSHOOTER data reduction pipeline version 
3.3.9] after which the spectrum was corrected for telluric ab- 
sorption using a synthetic transmission spectrum computed us- 
ing the molecfit packag¢?| fitted to the visual and near-IR data. 
The resulting spectrum is shown in Fig. [I] 


3. Stellar Parameters 


The XSHOOTER spectrum of HD 152384 (Fig.[I) is dominated 
by strong absorption lines from H1 and the Can H (396.8 nm) 


l http://www.eso.org/sci/software/pipelines/ 
2 http://www.eso.org/sci/software/pipelines/skytools/molecfit 
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Fig. 1. UV-Optical XSHOOTER spectrum of HD 152384 (black line) with the location of the most prominent emission lines indicated. Also 
shown (red line) is the spectrum of HD 207222 (Sp. Type = AO), obtained from the XSHOOTER Spectral Library (Gonneau et al.|2020), shifted 


for clarity. The inset shows the region around the Ha line (656.3 nm). 


Table 1. Stellar properties of HD 152384 derived in this paper. 


d Sp. Type Ter Vrad v sini Ly R, Mx Prot i Age 
[pc] [K] [km s7'] [km s7'] [Lo] [Ro] [Mo] [hours] P] [Myr] 
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log A [um] 
Fig. 2. Spectral Energy Distribution of HD 152384. Observed photom- 
etry (black dots) is compared to a{Kurucz| (1991) model for a AO (Ter 
= 9500 K) star (grey line), reddened using the |Cardelli et al.| (1989) 
normal extinction law and a colour excess E(B-V) = 0703. 


and K (393.4 nm) lines. From a Gaussian fit to these absorption 
lines we derive a heliocentric radial velocity for the star of —10 
+ 2 km s™!. We compare the spectrum of HD 152384 with a 
template spectrum for a AO star (HD 207222), taken from the 
XSHOOTER Spectral Library (Gonneau et al.!2020). The stellar 
absorption lines seen in this spectrum are a near-perfect match to 
the absorption lines seen in HD 152384 (Fig.|1), confirming its 
previous spectral classification of AO V ouch 1982 1982). We also 
derive a value for the stellar v sini of 116 + 20 kms™ from the 
FWHM of the Mg u 448.1 nm line (seen in absorption) using the 
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Fig. 3. TESS broad-band (600-1000 nm) normalized photometric time 
series of HD 152384 (black dots), covering the time period May 21, 
2019 to June 18, 2019. In addition to the periodic variations on a 8.286 
hour time scale, a flare can be seen in the data around 20.5 days. The 
red line shows the Gaussian Process fitted to the data. 


method of (2002), correcting for the lower spectral 


resolution of XSHOOTER compared to the instrument used by 
Royer et al. 


Using optical photometry retrieved from the SIMBAD 
database, as well as ultraviolet, near- and mid-IR data retrieved 


from the TD1 (Thompson et al.|/1978), 2MASS 


2003) and AIIWISE (Cutri et al.|/2013) catalogues, we con- 
struct the Spectral Energy Distribution of HD 152384 (Fig. [2). 


As is shown in this Figure, the optical and near-IR photome- 
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Fig. 4. Detected emission lines in the XSHOOTER spectrum of 
HD 152384. Velocity scales for all lines have been corrected for the 
radial velocity of —10 km s~! derived from the photospheric absorption 
lines. 


try of HD 152384 can be fit well by a mildly-extincted (Ay = 
009) model for the photosphere of a AO V star (Te = 9500 K), 
whereas a moderate (Ljp/L, = 0.05) excess emission above the 
expected photospheric levels can be seen longward of 2 microns. 
Using the integrated stellar flux from the fitted [Kurucz] (1991) 
model for the stellar photosphere, as well as a distance of 128.3 
+ 0.6 pe derived from the Gaia EDR3 parallax 


joration||2020), we derive a stellar radius of 1.9 + 0.2 Ro and 
stellar luminosity of 27.4 + 2.1 Lo for HD 152384. Using the so- 


lar abundance pre-main sequence evolutionary tracks of |Bressan| 
this combination of stellar temperature and lumi- 
nosity corresponds to a 2.1 + 0.1 Mo star located close to the 
Zero-Age Main Sequence, indicating an age larger than 4.5 x 
10° years. This lower limit to the stellar age of HD 152384 is 
compatible with membership of the Sco-OB2 association (e.g. 
2008). Our derived stellar parameters are 
summarized in Table 

An inspection of archival TESS data for 
HD 152384 show the star to be variable at the 0.1% level, with 
a strict periodicity (Fig. B). We estimate the period by fitting a 
Gaussian Process (GP) with a rotation kernel 
eq. 56), which takes into account the presence of cor- 
related noise in the TESS data. We find a period of 8.286 + 0.006 
hours, which we identify with the stellar rotation period. In ad- 
dition to the regular brightness variations which we attribute to 
stellar rotation moving a different pattern of starspots in and out 
of our line of sight, we also note that a small (0.2%) flare is 
present in the TESS data around 20.5 days, which could be due 
to further low-level stellar activity. Adopting a stellar radius of 
1.9 Ro, the rotation period found above corresponds to a rotation 
velocity of ~280 km s™!, or 62% of the stellar break-up velocity 
of 460 km s~!. Although relatively fast for a main-sequence star, 
this rotation velocity is in-line with those found for other young 


Table 2. Detected emission lines in the spectrum of HD 152384. 


Aab. [nm] Ident. EW [A] 
729.15 [Can] -0.20 
732.39 [Can] —-0.18 
849.81 Can —1.60 
854.21 Can -2.87 
866.21 Can —2.44 
838.78 Fer —0.14 
851.41 Fei —0.19 
882.41 Fei —0.10 
804.78  Mgı —0.15 
880.67 Mgr —0.14 
791.38 Sir —0.12 
807.46 Sii —0.09 
835.04 Sir —0.05 


early A-type stars (Wolff et al.|2004}|/Zorec & Royer|2012), pro- 
viding further support for the classification of HD 152384 as a 


young star. With the vsini value derived from the Mgu 448.1 
nm line, we also use the rotation velocity to derive a value for 
the stellar inclination of i = 24 + 4°. 


4. Circumstellar Material 


There are several indications for the presence of circumstellar 
material in HD 152384. The infrared excess emission seen long- 
ward of 2 um in the SED (Fig. B) suggests the presence of warm 
dust in the vicinity of the star. In addition, the XSHOOTER spec- 
trum displays a number of emission lines superimposed on the 
photospheric absorption spectrum. We identify these lines with 
emission from Can, Sir, Mgı and Fe1 (Table [2). All of the de- 
tected emission lines are centered on the stellar radial veloc- 
ity, and thus appear to be clearly associated with HD 152384. 
The emission lines display a clear double-peaked profile, with a 
peak-to-peak separation of around 70 km s~! (Fig. 4). 

The observed emission line profiles are reminiscent of the 
line profiles exhibited by the rotating gaseous disks seen around 
many young stars (e.g. |van den Ancker et al.[2004). These lines 
cannot be due to stellar activity, as (1) both the spectral type of 
AO V and the low level of variability seen in the TESS data sug- 
gest low levels of stellar activity, (2) stellar activity would result 
in emission components in the Cam H and K lines (e. g. [Vaughan] 
[et al.|[1978), which are not seen in our data and (3) stellar ac- 
tivity would not give the double-peaked emission line profiles 
seen in the XSHOOTER spectrum. The lines profiles are also 
unlikely to be solely due to a stellar wind as they do not exhibit a 
P Cygni like shape. We thus conclude that an origin in a circum- 
stellar disk is the only plausible explanation for the presence of 
the emission lines seen in the XSHOOTER spectrum. 

With the exception of the [Can] lines, which appear to be 
more symmetric, all detected emission lines display a distinct 
asymmetry, with the blue peak slightly higher than the red peak. 
This slight asymmetry could be a reflection of an asymmetric 
distribution of material within the disk, or be due to a wind from 
an optically thick disk (e.g. 
(2020). Alternatively, there could be some additional extinction 
in our line-of-sight towards the part of the disk which is mov- 
ing away from us. Higher spectral resolution data are needed to 
distinguish between these possibilities. 
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Following the approach in|Carmona et al.|(2011) we model 


the detected emission lines using a toy model for gas in a Kep- 
lerian orbit around HD 152384 assuming a flat disk, with inner 
radius Rin and outer radius Rou. The intensity of the emission 
decreases as J(R) œ (R/Rin)®, with R the radial distance from 
the star. The total intensity of the model is computed as a func- 
tion of velocity for a given inclination i and convolved with a 
Gaussian with width corresponding to the XSHOOTER spectral 
resolution in the visual arm. Best fit parameters were selected 
by comparing a grid of simulated spectra with an average of the 
continuum-subtracted Ca n emission lines using a y? minimisa- 
tion procedure. Examples of the resulting simulated spectra are 
shown in Fig. [5] 

Assuming that the disk has the same inclination as derived 
for the stellar photosphere (24 + 4°), the emission from the disk 
needs to originate close to the star (Rout < 0.3 au) to be able to 
reproduce the double-peaked emission line profiles seen in the 
spectra. However, we note that in our model the disk inclination 
is to some extent degenerate with the disk outer radius, so the 
emanating region would move correspondingly further out if the 
disk is seen under higher inclination than the value derived from 
the stellar photosphere. The lack of high-velocity wings in the 
observed line profile may indicate the presence of a small in- 
ner gap, which can be reproduced in our model by setting Rin 
to a value of around 0.03 au. Varying a, the exponent of the in- 
tensity distribution, has only a small effect on the width of the 
computed line profile. It is therefore poorly constrained by our 
observations. 

What makes the emission-line spectrum of HD 152384 
highly unusual is the complete absence of Hydrogen and He- 
lium emission (c.f. Fig. p. suggesting that the disk is highly de- 
pleted in volatile elements. To the best of our knowledge this 
spectrum with strong emission lines from refractory elements, 
but completely lacking emission from volatile species is unique 
within the group of young stars. However, similar spectra are 
seen in rare instances in a completely different class of stars: 


metal-polluted white dwarfs (Gänsicke et al./2006 


emission line systems also exhibit excess infrared emission in 
their SED caused by the presence of circumstellar dust. In the 
case of these white dwarfs the dust disk is the result of the tidal 


disruption of rocky planetesimals (Veras et al./2015}|/Malamud & 
Perets|2020}, and the double-peaked emission lines in their spec- 
trum is thought to originate from a gaseous component in these 


circumstellar debris disks. The exact mechanism generating the 
gas is still a matter of debate and the proposed theories include: 
dust sublimation at the inner edge of the debris disk followed 
by radial spreading E ond te 
crushing the debris into gas (Kenyon & Bromley|2017), and the 


presence of a dense planetesimal orbiting within the disc and 


disrupting the dust (Manser et al./2019). 


5. Discussion and Conclusions 


HD 152384 appears to be unique within the group of young 
intermediate-mass stars with disks. Its spectrum shows double- 
peaked Ca, Si, Mg and Fe emission lines, but no Hydrogen or 
Helium emission, indicative of a disk which, although rich in re- 
fractory elements, is highly depleted in volatile materials. The 
slightly higher blue peak seen in most lines could be due to a 
slight asymmetry in the disk, could indicate the presence of a 
disk wind, or that there could be some additional extinction to- 
wards the redshifted disk material. Assuming that star and disk 
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Fig. 5. Simulated emission lines from a flat circumstellar disk using 
the toy disk model. The grey lines show the continuum-subtracted aver- 
age Cau line detected in HD 152384, whereas the green lines show the 
model with the parameters i = 24° (corresponding to the stellar incli- 
nation), œ = 1.0, Ri, = 0.03 au, Rout = 0.15 au (corresponding to the 
best fit parameters). In each panel the red and blue lines illustrate the 
effect of varying one model parameter (top left: inclination, top right: 
a, bottom left: disk inner radius in au, bottom right: disk outer radius in 
au). 


are seen under the same inclination angle of ~ 24°, the width and 
peak-to-peak separation of the detected double-peaked emission 
lines would suggest an origin of these lines in a region close (ra- 
dius < 0.3 au) to the star. 

These properties of HD 152384, although unique for a young 
star, bear a striking resemblance to a rare subset of metal- 
polluted White Dwarfs, in which the presence of an emission- 
line spectrum composed solely of refractory elements is com- 
monly explained by the destruction of rocky planets 
[& Xul/2017). We speculate that similar processes could occur 
around a young star with a newly formed exoplanetary system, in 
which the disruption caused by planetary migration could cause 
collisions between rocky planets. We note that the compact (< 
0.3 au) disk around HD 152384 is smaller than the disks seen in 
most Herbig Ae/Be stars and references therein), 
but its location is similar to rocky planets in our own solar sys- 
tem. In addition, simulations of tidal disruption of giant plan- 
ets due to interactions with their host star show that remnants 
of exomoons can end up in orbits that are similar to the disk 
seen in HD 152384 (Martinez et al.|2019). Interestingly, the rel- 
atively low fractional infrared luminosity (Lip /L, = 0.05) places 
HD 152384 within the group of extreme debris disks defined re- 
cently by [Moér et al.| (2021), for which large amounts of warm 
dust may stem from recent giant impacts between planetary em- 
bryos during the final phases of terrestrial planet growth. Al- 
though cool molecular gas has been found at relatively large dis- 


tances in a number of these objects (Moor et al.2017 
2020), our analysis on HD 152384 constitutes the first de- 
tection of gas solely consisting of refractory elements around 


such an object, leading further support for the hypothesis of the 
destruction of planetary embryos. 

Signatures of selective accretion of circumstellar gas have 
also been found in the abundance patterns seen in the photo- 
spheres of A Boo stars and Herbig stars (Venn & Lambert|1990} 
2018), 
and may be due to Jupiter-like planets blocking the accretion of 
part of the dust. Several systems of similar age as HD 152384 
such as £ Pic also display variable absorption features in their 
spectrum which are interpreted as due to infalling evaporating 
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comets (e.g. 1987 1991) 
2014} | Vidal-Madyar et al./2017}|Rebollido et al.!2020), demon- 
strating the dynamic nature of newly formed planetary systems. 


Collisions are also believed to have played an important role in 
shaping our own solar system, and may have led to the formation 


of our Earth’s moon (e.g.|Canup & Asphaugh)2001). We postu- 
late that HD 152384 may currently be experiencing a similarly 


dynamic phase in its evolution. One prediction of our hypothesis 
of disruption by planetary migration is that HD 152384 should 
possess one or more giant planets at radii larger than a few tens 
of an au which would be responsible for the disruption of the 
inner rocky planets and the creation of the circumstellar disk 
reported here. High contrast imagers on the next generation of 
Extremely Large Telescopes will likely be able to detect these. 
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